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We investigate the effect of non-equilibrium and interactions on shot noise in v = 2 quantum 
Hall edges, where interactions between co-propagating edge modes are expected to give rise to 
charge fractionalization. Using the method of non-equilibrium bosonization, we find that even 
asymptotically the edge distribution function does not thermalize, but instead depends in a sensitive 
way on the interaction strength between edge modes. We compute shot noise and Fano factor from 
the asymptotic distribution function, and from comparison with a reference model of fractionalized 
excitations we find that the Fano factor can be close to the value of the fractionalized charge. 



In contrast to three spatial dimensions, where excita- 
tions of an interacting many particle systems often carry 
the same quantum numbers as in the non-interacting 
case, interactions in Id systems completely change the 
character of the excitation spectrum [1]. A prototype 
model for this physics is the Luttinger model, where 
electrons are no longer well defined quasi-particles, and 
where electronic excitations decompose into spin and 
charge parts moving with different velocities [2] 

An important example of interacting Id systems are 
the edge states of incompressible quantum Hall liquids 
[31 H] , where as a result of strong interactions charge frac- 
tionalization can occur [51-[TU] and manifest itself in shot 
noise pTI - fT4] . For the case of filling fraction v = 2, there 
are two co-propagating chiral edge channels. In the pres- 
ence of a short range interaction between them, a pulse of 
charge e injected into edge mode one at a first quantum 
point contact (QPCl) decomposes into a charge pulse 
and a neutral pulse. In the charge pulse, a charge e* 
travels on mode two and e/2 + \/ — 4(e*)2 on mode 
one [TS]. In the neutral pulse, there is a charge — e* on 
mode two and a charge e/2 — -^/e^ — 4(e*)^ on mode one. 
In this way, by exciting edge channel one via a partially 
transmitting QPCl, high frequency charge noise is gener- 
ated on edge mode two [5]. A second QPC2, allowing for 
partial transmission of channel two, transforms the high 
frequency charge noise into zero frequency shot noise [5] ■ 
If a fractional charge e* was scattered in the same way 
as an electron, one would expect a Fano factor e* for the 
charge noise generated at QPC2. 

There is an alternative way of looking at this prob- 
lem, related to the concept of energy relaxation [TBlI19j . 
Interactions play a crucial role in the thermalization pro- 
cess that drives a system through states described by the 
Gibbs equilibrium ensemble. Generically, the dynamics is 
only constrained by two integrals of motion, total energy 
and total particle number. Integrablc models like the 
V — 2 quantum Hall edge have infinitely many integrals 
of motion, and therefore it is not clear if an equilibrium 
state can ever be reached [20]. If the two edge modes 
are driven out of equilibrium with respect to one an- 
other, the system relaxes towards a non-thermal steady 




FIG. 1: (color online) Sketch of a = 2 Hall bar with a 
QPCl, where inner modes ("2" , light blue) are fully reflected, 
while partial transmission of outer edge modes (" 1" , black) 
is possible. At QPC2, the opposite situation is realized. The 
shaded area is the interaction region, where partial energy 
relaxation takes place. The upper ede is biased with voltage 
V at contact 1, current noise is measured at contact 3. 



state, whose distribution function determines shot noise 
at QPC2. The corresponding Fano factor depends on the 
strength of the interaction between the two edge modes, 
and in general neither agrees with the fractional charge 
e* introduced above, nor with the result for two equili- 
brated edge modes. For the special case of a half open 
QPCl however, the Fano factor is close e*/e, suggesting 
an interpretation in terms of charge fractionalization. 

We consider the setup Fig. [l] where a Hall bar is 
pinched by two QPCs. The outer edge mode is labeled 
"1" and the inner one "2". The top and bottom edges 
originate at zero temperature from reservoirs at voltages 
Vi = V and V2 = 0. At QPCl, the outer modes are 
partially transmitted with probability a, while the in- 
ner ones are fully transmitted; as a consequence, only 
the outer mode becomes noisy. After QPCl, the two 
edge modes interact over some distance (shaded area in 
Fig. [1] with length £) before reaching QPC2. Here, the 
outer modes are fully reflected while the inner ones are 
partially transmitted with probability p. Current noise 
is then measured at contact 3. Using the recently de- 
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veloped non-equilibrium bosonization technique [2T1 [55] 
within a quantum-quench model [23H25] , we compute the 
shot noise at QPC2, with particular emphasis on its de- 
pendence on the strength of the interaction between the 
two edge modes. 

The edges and QPC2 are described by the following 
Hamiltonian ( if not otherwise stated, we adopt units 
with h = ks — i ) '■ 

Jfj, =271 [vipj^^{x) + V2pI,jj{x) + Vi2Pl,r,{x)p2,ri{x)) 

J X 

J^QPC2 = t2lpl^u{x)^2A^) + ^-C- (1) 

Here, describes chiral modes with velocities ui, W2, 
and with a local interaction vi2 between them. J^qpc2 
describes scattering of electrons at QPC2 with t2 the 
scattering amphtude. The fields Pi,r,{x) in ([T]) describe 
density fluctuations, they are related to bosonic dis- 
placement fields by pi{x) = dx(f>i{x) /2tt; here "i" la- 
bels different edge modes, and -q = u,d the upper 
and lower edge. The bosonic fields have commuta- 
tors [(l){x) , (j){y)] = Z7rsign(x — y), and the fermionic 
field is represented as 4'riix) = (27rQ;)~^/^ e*'^'''^^^ with 
a denoting a short distance cutoff on the scale of the 
magnetic length. For later reference, we decompose 
the bosonic fields as (t>n{x) = frjix) + ^pljix), f,-i{x) = 
E,>ov/W9^e-9"/2e*-''9-6^(s^g), where = ±f re- 
spectively for right (u) and left (d) movers, and h\h are 
canonical bosonic operators. 

Following [TH], we do not model QPCl explicitly but 
instead consider its effects on the downstream electron 
distribution of mode (Iw), which is modeled as a "double 
step" function 

/(e) = a e{-e + ^i) + (f - a) e{-e + /i2), (2) 

where /ii = {l—a)eV and ^2 — —a eV (assuming eV > 0) 
are chosen such that the average density in mode (lu) 
corresponds to zero bias. This choice for / is made to 
ensure that the interaction V12 does not give rise to a 
density shift in mode (2m). 

Next, we consider the effects of the inter-mode inter- 
action on the distribution functions by formally turning 
on the interaction right after QPCf . We use the model 
of a quantum quench, where the interaction is suddenly 
turned on at time t = 0, thus changing the coupling 
constant V12 from zero to a finite value. The full inter- 
acting problem can now be diagonalized by means of a 
Bogoliubov transformation M. For co-propagating states 
{viV2 > 0), M can be represented by the following rota- 
tion matrix in operator space : 

M=( ^"^M, (3) 

\^ — sm 6^ cos & J ^ ' 

allowing to express Jif in terms of new (rotated) fields 
A.g = J2j Mijbj_q. The mixing angle 6 expresses the 
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FIG. 2: (color online) Steady state distribution of edge mode 
{2u) asymptotically far away from the QPCl. (black full line) 
Nonequilibrium distribution obtained from Eqs. ( |7|9[ ) by sum- 
ming over all cumulants. (blue dash-dotted line) distribution 
obtained by retaining only the gaussian term, (red dashed 
line) fully equilibrated distribution at effective temperature 
T* = eVy/{3/2)a(l - a)/-K. The mixing angle is 61 = 0.47, 
and the transmission probability of QPCl is a = 1/2. 

strength of the interaction through the relation tan 20 = 
Vi2/{vi — V2)- At this point the rotated operators evolve 
in the Heisenberg picture as Piq{t) = e~*''"'*/3iq(t = 0), 
where the new velocities are given by Wi(2) = "^^1(2) cos^ 0+ 
V2(i) sin^ 6 ± ^vi2 sin 29. 

As a final step, we undo the Bogoliubov transformation 
in order to express the l3iq{tQ) in terms of the original ba- 
sis. As a result, we obtain a relation between the bosonic 
operators at > and those at t = : 

blq{to) = Uq{to)blq + Sq(to)&2g (4) 
hqih) ^ Sq{to)blq + Wq(io)&2g 

where biq = biq{t = 0). Now all the time dependence is 
encoded in the coefficients 

Uq{to) = cos^ 9 e-'"^^'" + sin^ e e'"^^''° (5) 
Vq{to) = cos^ 9 e-'"^^'" + sin^ 9 e-"^^''" 

where "fg = sin 29. To leading order in the tunneling 
amplitude t2, the noise in the tunnel current at QPC2 can 
be expressed in terms of greater (lesser) Green functions 

G>^<^(e) m as 

= ^^^G< (6)G>(e) + G<(6)G>„(e),(6) 

with G^(e) = G^(— e). Using the boson representation 
of electron operators, we can compute G^J*^'' (r) of the 
fully interacting edge mode. Due to the non-equilibrium 
distribution of edge mode (lu), calculating the expec- 
tation value of a product of bosonic exponents is highly 
non trivial. Here we discuss the results for the "long time 
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FIG. 3: (color online) Shot noise after QPC2 as a function of 
transparency a of QPCl, normalized to one at a = 1/2, for 
a mixing angle 6 — 0.47. Comparison of full non-equilibrium 
result (full black line) , reference noise of noninteracting elec- 
trons (dashed red line) , and noise in a fully equilibrated ther- 
mal state (dash-dotted blue line). 



limit" of the Green function, in which the system reaches 
a non-equilibrium steady state: 

G2uir) = {^Lito + T,x)4>2u{to,x)) (7) 

1 1 1 



Go<(r) 



Here, Gq (r) is the equilibrium Green function of edge 
mode 2 in the presence of interactions. All the infor- 
mation about non-equilibrium effects is contained in the 
average over bosonic coherent states in Eq. ([7|, where 
A(g,io,T) = iSq(io)(27r/gL)i/2e ' 
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As emphasized in [21] . non-equilibrium effects make 
the theory non-Gaussian, and higher order cumulants ap- 
pear in the evaluation of the above expectation value. In 
order to compute the expectation value over the non- 
equilibrium state, we refermionize the bosonic operators 
introducing new fermionic operators |27| : 



6L ^ (27rM)V2^c1„_,+,ci„„ 

k 



(8) 



Since the bosonic operators describe free particle-hole 
excitations, also the c-operators are free and therefore 
can be connected to the incoming states via a scattering 
matrix. Then, the desired expectation values of products 
of Fermi operators can be evaluated by using Wick's theo- 
rem and an appropriate fermionic density matrix pi„.The 
crucial step now consists in noticing that the computa- 
tion of higher order cumulants is similar to the problem of 
full counting statistics, and using Klich's trace formula 
[551 HH] they can be expressed in terms of a Fredholm 
determinant of the Toeplitz type, normalized to its zero 



FIG. 4: Fano factor as a function of transparency of QPCl 
for the mixing angle 6 = OAT. At a = 1/2 the Fano factor is 
F = 0.45. 



temperature, equilibrium value 

dot [1 + (e-'"^- 



(9) 



dot [1 + {e-'^- - l)6l(-e)] ' 

where /(e) is given by Eq. ([2|. The scattering phase 
dr = TTjgujr(to^O) Contains information about the inter- 
edge interaction, and the window function 

UJr{tQ,0) = 9{x - {Vi + V2)tQ) - 9{X - {Vi + V2)tQ - V2t) 

- 9{x ~ 2i2to) + 9ix - i2i2to + t)) (10) 

restricts the action of the determinant over a time range 
T. The window function represents a combination of 
two square pulses which separate on a time scale t <C 
£/{v2 + vi). In this limit, the expectation value of bosonic 
coherent states factorizes into a product of two deter- 
minants having the same scattering phase, and we can 
rewrite Eq. ^ as G^„(t) = G^{t)AI{S). The determi- 
nant Eq. ^ can be evaluated numerically by carefully 
defining a regularization scheme j22j . 



Finally, the lesser Green function G^^(e) 

/,--,7fl/2-.l-7fl/2 



9{—e)/vi V2 is easily evaluated due to its equi- 

librium nature. Fourier transforming Eq. ([7| into en- 
ergy space, we can compute the distribution function 
at QPC2; as a consequence of interactions, the distribu- 
tion function is broadened from a single step (see Fig. [2| . 
However, it does not have the same functional form as 
a Fermi distribution, but rather describes a true non- 
equilibrium steady state. The distribution obtained by 
only retaining the Gaussian term in the cumulant ex- 
pansion clearly deviates from the full one, making evi- 
dent the necessity for including higher order terms. The 
non-equilibrium distribution also deviates from an equi- 
librium Fermi distribution with effective temperature 
T* = eVy/{3/2)a{l - a)/7r, obtained by assuming that 
the two edge modes fully equilibrate and that each of 
them carries half the energy flux injected into the upper 
edge via QPCl. 

In Fig. (|3|, we display the dependence of low frequency 
noise computed from Eq. (|6| on the transmission a of 



4 




FIG. 5: (color online) Upper panel: Fano factor as a func- 
tion of the mixing angle for transmission a = 1/2 of QPCl. 
(red, dashed) Fully equilibrated edge, F is independent of 
interactions, (black, dotted) Full non-equilibrium situation, 
(blue, continuous) Reference model diluted system of frac- 
tional charges {F = (l/2)sin20). Lower panel: Fitting the 
non-equilibrium noise displayed in Fig. 3 by a function propor- 
tional to (a(l — a)Y , the parameter d depends on the mixing 
angle 9. 



QPCl, normailizing the noise by its value at a = 1/2. 
One clearly sees that it deviates from both the standard 
free fermion dependence a(l — a), and from the equilib- 
rium result with ^J7L(l^^~a). 

To put the strength of the noise at QPC2 in per- 
spective, we define a reference noise expected for non- 
interacting electrons tunneling through both QPCl and 
QPC2 along a single edge, obtained by using the distri- 
bution Eq. (2) in Eq. (6) 



S,ci{ijj^{))=Aepa{\-a)I with I = ■ (H) 

Here, the transmission probability of QPC2 is p = 

7.= 



\t2 



Defining a Fano factor F = S/S^cf as the ratio be- 
tween the noise obtained from Eq. ^ and the reference 
noise Eq. ( [Tl] ) (see Fig.|4|, we can make contact with the 
concept of fractional charges described in the introduc- 
tion. Assuming that the non-equilibrium state of mode 
(2m) could be described by a diluted system of fractional 
charges e* = (e/2) sin 20 15 , it is plausible that the noise 
would be given by Eq. (11 1 when replacing the charge e 



by e* , giving rise to a Fano factor e* /e (strictly speaking, 
there are two charges e* and — e* associated with each 
electron, but on the other hand the " currentless" noise of 



factor is shown as a function of mixing angle for the spe- 
cific transmission a = 1/2 of QPCl. For this value of a, 
there is a surprisingly good agreement between the value 
e*/e— (1/2) sin 20 [TF and F of the fuU non-equilibrium 
noise, suggesting that the Fano factor can be interpreted 
as being due to formation of fractionalized charges in the 
1^ — 2 quantum Hall edge. 

We find that the zero frequency noise power depends 
in a singular way on a in the limit a <C 1 |29| . To ob- 
tain the noise in this limit, the functional determinant 
can be approximated by its long time asymptotics (valid 
for eVr > 1) Ar{S) ~ exp(-|r|/(2r^)), where the de- 
phasing rater^^ = -(ey/27r) log[l-4a(l-a) sin^(5/2)]. 
Knowledge of Ar{S) for large times allows to accurately 
calculate the distribution function of mode (2u) for en- 
ergies e ^ eV. However, for a ^ 1 the distribution 
function only deviates from a step function on the scale 
aeV, such that the long time asymptotics allows an exact 
calculation of the distribution function. Using Eq. ^ 
and taking the a ^ 1 limit, we find S{lu — > 0) ~ 
8palog(l/a) sin2((5/2)ey(eV/i7r2) . 

A useful way to characterize the nonlinear dependence 
of experimentally measured shot noise on the transmis- 
sion probability a of QPCl is by fitting it to a function 
proportional to (a(l — a))'' [301 • For the reference noise 
of Eq. (11), d is trivially equal to unity. For "thermal" 



noise generated by the effective equilibrium temperature 
T* , one finds d = 0.5. For the full equilibrium noise, we 
find that its dependence on a can be well fitted by the 
above power law, and that d varies from d = 0.85 for 
e = 7r/16 to d = 0.68 for 9 = 7r/4, see Fig.jsj In this way, 
from knowledge of d the mixing angle can be inferred, 
without making use of the Fano factor. 

In summary, due to the joint effect of interactions and 
non-equilibrium, the distribution function of an originally 
unbiased, zero temperature mode (2u) interacting with 
a noisy mode (lu) evolves towards a non-thermal steady 
state that depends on the interaction strength in an es- 
sential way. Comparing the shot noise and Fano fac- 
tor from our numerically exact calculation with a simple 
model of charge fractionalization, we find that the Fano 
factor can indeed be interpreted in terms of charge frac- 
tionalization in the v = 2 quantum Hall edge. 
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Eq. ( 11 ) already contains contributions of both electrons 
and holes, so at least for the specific value of a = 1/2 this 
would amount to double counting) . In Fig. [s] , the Fano 
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